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EFFECT OF ELfCTROTHERMAL INSTABILITIES ON BRAYTON- AND RANKINE- 
CYCLE MAGNETOHYDRODYNAMIC SPACE-POWER GENERATION SYSTEMS 
by Allan R. Bishop and Lester D. Nichols 
Lewis Research Center 

SUMMARY 

The effect of fluctuations in the electron density on Brayton- and Rankine- cycle 
magnetohydrodynamic generators is analyzed. The generators are designed for mini- 
mum magnet volume, the Bra 5 d;on cycle using neon seeded with cesium and the Rankine 
cycle using lithium seeded with cesium. 

An important consequence of the electron density fluctuations is an increase in Joule 
heating within the generator and a reduction in output power. The power loss increases 
as the maximum cycle temperature decreases. An increase in the entrance magnetic 
field can compensate for the reduction in output power. For a given entrance magnetic 
field, generator length, generator efficiency, and maximum cycle temperature, the 
power output of the Brayton- cycle generator is maximum at an entrance Mach number 
of 1.0. 


INTRODUCTION 

The need for increasingly larger, lighter sources of power in spacecraft has stim- 
ulated interest in alternatives to the conventional turbogenerator set. A closed- cycle 
magnetohydrodynamic (MHD) generator using h^h-temperature gases is one possibility, 
and has been investigated by several authors (see refs. 1 and 2). Nichols (ref. 3) com- 
pared two different MHD cycles with a turbogenerator set and studied the range of pa- 
rameters which give comparable system weights. 

All of these references use the conventional formulation for electrical conductivity 
and Hall parameter in the ionized gas. This formulation neglects any nonuntformities in 
the plasma. However, some investigators (refs. 4 and 5) have noted that in the presence 
of strong magnetic fields nonuniformities develop in the plasma which lead to a reduction 
in the electrical conductivity and changes in other properties. Analytical studies of this 



phenomenon have developed the concept of a critical Hall parameter. When the micro- 
scopic Hall parameter is larger than this critical value, any variation is amplified and 
appreciable fluctuations in the local properties occur. The value of this critical Hall 
parameter depends on the electron temperature and, to a lesser extent, on the amount of 
radiation and thermal diffusion present. Solbes (ref. 4) has derived anal 3 rtical expres- 
sions for the critical Hall parameter, the changes in bulk properties, and the magnitude 
of the fluctuations. 

In this report, the ej^ressions developed by Solbes (ref. 4) are used to modify the 
calculations given by Nichols (ref. 3). The same thermodynamic cycles discussed by 
Nichols are considered, and similar assumptions are made. The only change is in the 
formulation of electrical conductivity. The effective value predicted by Solbes is sub- 
stituted whenever the Hall parameter is above its critical value. This change should 
make the estimates of generator performance more accurate and the comparison between 
generators more valid. 


SYMBOLS 


B magnetic field strength 
eVj^ ionization potential of the seed 
j current density 

K parameter accounting for variations in electron temperature and loss factor 

k Boltzmann constant 

M entrance Mach number 

m^ atomic mass of the carrier 

m^ electron mass 

mg atomic mass of the seed 

in average particle mass 

N electron number density 

N total gas (seed plus carrier) number density 

o 

N total seed (ion plus neutrals) number density 
s 

electron-neutral collision cross section of carrier 

Q_ electron-neutral collision cross section of seed 
s 

r parameter specifying importance of Coulomb collisions 
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S seed fraction 

T bulk gas temperature 

Tg electron temperature 

Z parameter related to effective conductivity 

a degree of ionization, N^/Ng 

^ local Hall parameter 

^crit critical Hall parameter 
^eff effective Hall parameter 

A composite loss factor 

V electron collision frequency, v . + i/ 

electron- ion collision frequency 
electron- neutral collision frequency 
cr local conductivity 

CTeff effective conductivity 

( ) average of quantity over space 


GENERATOR CYCLE SPECIFICATION 

The following discussion is based on the report by Nichols (ref. 3), and a more de- 
tailed analysis can be found there. Two magnetohydrodynamic generators were consid- 
ered; a Rankine cycle using lithium seeded with cesium and a Bra 3 d;on cycle usii^ neon 
seeded with cesium. A turbogenerator using potassium was the reference cycle. Since 
waste heat must be radiated away in space, and the radiator is one of the more massive 
components as well as one that is sensitive to assumed operatir^ conditions, the two 
MHD generators were restricted so that all three cycles require equal radiator area. 

The reference cycle is a Rankine cycle using potassium and a turboalternator. A de- 
tailed description of the reference cycle is given in reference 3. 

The two MHD generators were assumed to be ideal except for Joule-heatii^ losses. 
No heat diffusion or radiation losses were considered. Both generators had a fixed 
length and a fixed aspect ratio. The values of j /o and jB were required to be con- 
stant over the length of the generator. This last restriction provides for uniform dissi- 
pation and uniform power generation along the duct. The generator area variation which 
minimizes generator volume at constant pressure ratio was assumed. The Brayton cycle 
had a pressure ratio of 0. 5, the Rankine cycle 0. 1. The efficiency of this minimum- 
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volume generator is very close to that of the maximum- efficiency case computed by 
Pleshanov (ref. 6). Following the procedure of reference 3, the maximum pressure, 
entrance Mach number, entrance current-density ratio, entrance magnetic field, and 
entrance Hall parameter for each cycle were found as a function of the maximum tem- 
perature. The requirement of equal radiator areas in both the MHD generators and the 
reference cycle specifies the generator efficiency as a ftmction of maximum cycle tem- 
perature (fig. 10 in ref. 3). From this generator efficiency the values of entrance Mach 
number and entrance current-density ratio were calculated (fig. 12 in ref. 3), based on 
the minimum- volume solution. The entrance magnetic field necessary to give the proper 
pressure drop within a generator length of 0. 5 meter was then determined. The gener- 
ator efficiency is the actual change in enthalpy divided by the isentropic chaise between 
the same pressure limits. 


CALCULATION OF PLASMA PROPERTIES 

The temperature of the electrons in the plasma may be h^her than the bulk temper- 
ature of the gas because of Joule heating. In steady state, the energy gained by the 
electrons through Joule heating must be lost to the heavy species by collisions. The 
balance can be written as 


^2 m„ 

i_= 3^1/ N k(T - T)A (1) 

a m 


This equation holds locally in the plasma and the value of each variable, particularly 
Ng, may fluctuate. 

In order to find an effective value for a, equation (1) is averaged over space 


< 


,2 


3 ^ k(T 
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m 

= 3 ^ k ( i'gN (T 
m e e e 


T)A) 


and an effective a is defined as 
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a =il>i 

eff / o\ 
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Then substituting for ( j /a) results in 




^eff 


m 




( 2 ) 


o 

The requirement that jB and j /a be constant is replaced by the requirement that (j> B 
2 

ajid (j) /cTgff constant aloi^ the generator. 


ASSUMPTIONS ABOUT FLUCTUATIONS AND PLASMA CONDITIONS 

The analysis of Solbes (ref. 4) is used to evaluate and to perform the averages 

in equation (2). An analysis similar to that of Solbes (ref. 4) has been made by Louis 
(ref. 5). However, Louis’s results are restricted to plasmas where Coulomb-type colli- 
sions are dominant. Solbes* s ejq)ression is more general in this respect, the relative 
influence of Coulomb collisions beii^ arbitrary. For this reason, Solbes’s development 
is followed rather than Louis’s. 

The fluctuations are assumed to be plane waves whose wavelength is much smaller 
than the dimensions of the channel. This means that the minimum dimension of the gen- 
erator must be about 10 centimeters (ref. 7). The bulk properties of the plasma, such 
as velocity, pressure, and temperature, are assumed to be uniform on a length scale 
which is large when compared with the wavelei^h of the disturbance. In the averaging 
process these properties are assumed to be constant, and only fluctuations with small 
wavelengths satisfying this requirement are considered. The fluctuations are also as- 
sumed to be periodic in space and to have small amplitudes. Only first-order terms in 
density fluctuations are retained, with one exception. Second- order effects must be 
considered in the Joule-heatii^ term of the energy equation since it is the source of the 
instability. The plasma consists of a weakly ionized seeded gas. 


Determination of Critical Hall Parameter 

An equation for the energy in the fluctuations is found by averaging the local energy 
equation and subtracting this averaged equation from the local energy equation. The 
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point at which the energy in the fluctuations changes from a decreasing function of time 
to an increasing one is the point where instability occurs and the fluctuations may be 
amplified. This point is reached when the Hall parameter is equal to, or greater than, 
the critical Hall parameter The value of as derived by Solbes is given by 

= 2V(1 + K)(r +K) (3) 

The parameter r is a measure of the importance of Coulomb-type collisions and varies 
between 0 and 1. Its value is given by 


r = 


ei 


V + V ■ 

en ei 



The parameter K relates the variation in electron temperature T^ and in loss factor 
A to the fluctuation in electron density. Its value is given in reference 4 by 


where 


K = 1 2 - c 

2 1 - a\eVj 





2 /m ' 

r(l + y) + (1 - r) ( — + y 

\“c / 


Nsd - Qs 

Ng-N3 Qe 


As shown by Solbes, the fluctuations in T^ are much smaller than those in un- 
less the degree of ionization a approaches 1. In this report, a is typically less than 
10 . Any fluctuations in T are therefore neglected. 


Calculation of Effective Conductivity and Hall Parameter 

Once the magnitude of is found, the effective conductivity is determined 

(ref. 4, eq. (3)): 
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Oeff = <a> 


( 4 ) 


2 2 ”il/2 

<1 - 

L(1 - r)2+ <^>2j 


^ (g) 

1 + z 


where (a) and (/3) are average values of the conductivity and the Hall parameter. Equa- 
tion (4) defines a new parameter Z which relates and (a> . 

The value of ^ used in equation (4) is the usual value based on the strei^h of the 
magnetic field and electron temperature. However, in experiments where density fluc- 
tuations occur, a different "effective" value of )3 is measured. A change in /3 is also 
indicated by the analysis. The effective Hall parameter is the open- circuit voltage 

perpendicular to the average current divided by the voltage parallel to the average cur- 
rent. Its value is given by 


^eff “ 


( 0 > 

1 + Z 


— ) 
1 + zy 


1 - r 

< 0 ) 


+ 



( 5 ) 


Calculation of Electron Temperature 


Now the value of T^ can be found by substituting equation (4) into equation (2). 


T 


e 


T + 



(1 + Z)m 

A) 


(6) 


where 



m 


e 



T 


e 


T, (j)^(l + Z)m 
3ke2(j.^Ng A> 



( 7 ) 


Note that variations in T^ have been neglected, so that the term T^ - T is a constant 
during the averaging process. 

Equation (7) is analogous to equation (39) in reference 3. Once the electron temper- 
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ature is determined, the calculation of the entrance magnetic field strength may be 
completed. 


Seed Fraction Specification 


The seed fraction (ratio of seed density to carrier density) remains to be specified. 
Consistent with the previous effort to minimize the weight of the magnet, the seed frac- 
tion which minimizes the entrance magnetic field is selected. In reference 3, a differ- 
ent criterion is used to specify the seed fraction, but the values used in both that work 
and this one are similar. 




(b) Rankine cycle boiling at 1430 K (2. TxlO'* N/ m^); working fluid, lithium seeded 
with cesium. 


Figure L ** Magnetic field strength required to provide a generator length of 0.5 meter 
as a function of seed fraction with entrance Mach number as a parameter. Maxi- 
mum cycle temperature, 2100 K. 
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Li figure 1(a) the entrance magnetic field of the Brayton cycle is shown as a function 
of seed fraction with entrance Mach number as a parameter. The maximum tempera- 
ture, maximum pressure, and current-density ratio are fixed. Note that the value of 
seed fraction which minimizes magnetic field is dependent on Mach number. At 
M = 0.5, it is about 0.006; at M = 1.0, about 0.004; and at M = 2.0, about 0.003. hi 
reference 3 a similar seed fraction of 0. 001 is recommended. 

For the Rankine cycle, figure 1(b) shows the entrance magnetic field as a function of 
seed fraction with entrance Mach number as a parameter. The maximum temperature, 
maximum pressure, and entrance current-density ratio are fixed. For M = 2. 0, the 
minimum magnetic field occurs at a seed fraction of 0.006; for M = 1.0, at 0.015. The 
curve for M = 0. 5 is very flat, giving a minimum value of magnetic field at a seed frac- 
tion of about 0. 050. However, increasii^ the seed fraction from 0. 01 to 0. 05 decreases 
the magnetic field by only 10 percent. A separate calculation shows that a decrease in 
output power of about 10 percent also accompanies this increase in seed fraction. The 
increase in seed fraction lowers the speed of sound in the plasma and reduces the input 
energy at constant entrance Mach number. A seed fraction of 0. 01 is used in the 
Rankine cycle operating at an entrance Mach number of 0. 5. 


RESULTS AND DISCUSSION 
Brayton Cycle 

The results of the Brayton- cycle calculations are shown in figure 2. The generator 
ler^h is 0. 5 meter, a specific radiator area equal to that of the reference cycle is 
maintained, and the optimum seed fraction is used. The entrance magnetic field 
strength is plotted as a function of maximum temperature, with Hall parameter, effec- 
tive Hall parameter, and maximum pressure as parameters, hi f^ure 2(a) the entrance 
Mach number is 0. 5, and the seed fraction is 0. 006; in figure 2(b) the entrance Mach 
number is 1. 0, and the seed fraction is 0.004; and in figure 2(c) the entrance Mach 
number is 2.0, and the seed fraction is 0.003. Figure 2 is similar to figure 20 in refer- 
ence 3, except that curves of have been added. 

The variation of the parameters in figure 2 is qualitatively similar to those in ref- 
erence 3. The entrance magnetic field increases with increasing pressure, increasing 
Hall parameter, decreasitg temperature, and decreasiig Mach number. Note the 
changes in with entrance Mach number. At high Mach number, ( 3 ^^^ is only a 

slight function of maximum temperature, while at low Mach number is a very 

strong fimction of temperature. As the maximum cycle temperature increases, 
approaches /8. At the point where )3 == the fluctuations are absent. For 
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Entrance magnetic field strength, B, tesla 


O 



Maximum cycle temperature, T, K 


(a) Entrance Mach number, 0.5; seed fraction, 0.006. (b) Entrance Mach number, 10; seed fraction, 0.004. (c) Entrance Mach number, 2.0; seed fraction, 0.003. 

Figure 2. -Entrance magnetic field strength required to provide a generator length of 0.5 meter as a function of maximum cycle temperature with maximum cycle pressure as a parameter, Brayton cycl^ working fluid neon seeded 
with cesium. ' 



^ ^ ^crit’ ^eff " ^ ^eff “ occurs at about 2400 K for M = 0. 5, as shoAvn in 

figure 2(a); and at temperatures above 2500 K for the higher entrance Mach numbers. 

The effect of instabilities is to increase the necessary entrance magnetic field, par- 
ticularly at the lower temperatures. This is shown in figure 3, where a comparison is 
made with the calculations given in reference 3. The maximum pressure is fixed at 
8x10^ newtons per square meter. The solid curves are taken from fi^re 2, with opti- 
mum seed fraction, while the dashed curves are from a calculation similar to that given 
in reference 3 with a seed fraction of 0. 001. 



Figure 3. - Effect of instability on the required magnetic 
field as a function of maximum cycle temperature with 
entrance Mach number as a parameter. Brayton cycle; 
maximum cycle temperature, 2100 K; maximum cycle 
pressure, 8x10^ newtons per square meter; working 
fluid, neon seeded with cesium. 
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Rankine Cycle 


Li figure 4 the results of the Rankine-cycle calculations are given. A generator 
length of 0. 5 meter, equal specific radiator areas, and optimum seed fraction are main- 
tained. Two boiling temperatures for the lithium are presented (1430 and 1645 K), and 
entrance Mach number is plotted as a parameter. The entrance magnetic field is a 
strong function of temperature, increasing rapidly as the temperature falls. The vari- 
ation with entrance Mach number is not as definite, depending on both boiling tempera- 
ture and maximum temperature. v 



Figure 4. - Entrance magnetic field strength 
required to provide a generator length of 
0. 5 meter as a function of maximum cycle 
temperature with entrance Mach number as 
a parameter. Rankine cycle,- working fluid, 
lithium seeded with cesium. 
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The comparison with reference 3 is shown in figure 5 for boUing at 1430 K. Again, 
the instabilities increase the necessary entrance magnetic field at a given maximum 
temperature. The solid curves are from f^ure 4, while the dashed curves are from 
reference 3 with a seed fraction of 0. 01. 

It is clear that for a fixed maximum cycle temperature, the effect of the electron 
density fluctuations is to increase the necessary entrance magnetic field. Alternatively, 
for a fixed entrance magnetic field, the maximum temperature in the cycle must be 
higher in order to produce the equivalent power. Both the Brayton cycle and the Rankine 
cycle show this increase. 



Maximum cycle temperature, T, K 


Figure 5. - Effect of instability on the required magnetic field as a 
function of maximum cycle temperature with entrance Mach num- 
ber as a parameter. Rankine cycle boiling at 1430 K (2. 7x10^ N/ m^); 
working fluid, lithium seeded with cesium. 
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Comparison of Cycles 


In figure 6 the output power at a fixed entrance magnetic field of 10 tesla is plotted 
against maximum cycle temperature with entrance Mach number as a parameter. Simi- 
lar data from reference 3 are also given. The higher cycle temperature necessary to 
maintain the same power output is evident in this graph. The maximum power is pro- 
duced by the Brayton cycle operating at M = 1. 0. Slightly less power is produced at 
M = 0. 5 and at M = 2. 0. The output power of the Rankine cycle is considerably lower 
than that of the Brayton cycle. The two boiling pressures selected for the Rankine cycle 
appear as discrete points, while the continuous range of maximum pressure in the 
Bra 3 d;on cycle forms a smooth curve. 

All of these results are qualitatively similar to those in reference 3. The fluctua- 



Figure 6. - Generator output as a function of maximum 
cycle temperature with entrance Mach number as a 
parameter. Entrance magnetic field strength, 10 tesla; 
Brayton cycle working fluid, neon seeded with cesium; 
Ranking cycle working fluid, lithium seeded with cesium. 
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tions in electron density change the numerical value of the results, primarily by de- 
creasing the conductivity and increasing the power dissipated in Joule heating. This im- 
plies either a lower power output or an increase in magnetic field strength to compensate 
for the loss. 


CONCLUSIONS 

The effects of electron density fluctuations on the performance of Brayton and 
Rankine magnetohydrodynamic (MHD) generators has been studied. The analysis indi- 
cates that 

1. The electron density fluctuations decrease the power output for a given generator 
and magnetic field strength. The power loss increases as the maximum cycle tempera- 
ture decreases. For the Brayton cycle, maximum power output is achieved at a Mach 
number of 1. 

2. Increasing the magnetic field strength will compensate for the output power loss 
due to fluctuations. Although a somewhat larger magnet is necessary, this does not 
limit MHD generators for space-power applications because the magnet weight is ex- 
pected to be a small fraction of the total system weight. 

3. There is a seed fraction which minimizes the entrance magnetic field strength 
for each generator. However, in some cases the advantage of a decrease in magnetic 
field is offset by a decrease in power output. 

4. The effective Hall parameter is considerably lower than the Hall parameter 
based on magnetic field and electron temperature. This means that the Hall voltage be- 
tween segments of the electrode is smaller, 

5. The power output of the Brayton cycle is larger than the Rankine cycle for the 
same maximum temperature, radiator area, length, entrance Mach number, and mag- 
netic field strength. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 7, 1969, 

129-02. 
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